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Q ABSTRACT 

Context. Using new homogeneous Luminosity Functions (LFs) in the FUV (VVDS) and in the FIR Hersciiel/PEP and 
Herscfiel/HerMES, we study the evolution of the dust attenuation with redshift. With this information in hand, we are able to es- 
timate the redshift evolution of the total (FUV + FIR) star formation rate density (SFRDtot)- By integrating SFRDtot, we follow the 
mass building and analyze the redshift evolution of the stellar mass density (SMD). 

Aims. This letter aims at providing a complete view of star formation from the local universe to z ~ 4 and, using assumptions on 
earlier star formation history, compares this evolution to what was known before in an attempt to draw a homogeneous picture of the 
^^ global evolution of star formation in galaxies. 

Methods. 

Results. The main conclusions of this letter are: 1) the dust attenuation Apuv is found to increase from z = to z ~ 1.2 and then 

starts to decrease up to our last data point at z = 3.6; 2) the estimated SFRD confirms published results up to z ~ 2. At z > 2, we 

■^l- observe either a plateau or a small increase up to z ~ 3 and then a likely decrease up to z = 3.6; 3) the peak of Apuv is delayed with 

f— ^ respect to the plateau of SFRD TO t and a likely origin might be found in the evolution of the bright ends of the FUV and FIR LFs; 4) 

using assumptions (namely exponential rise and linear rise with time) for the evolution of the star formation density from z = 3.6 to 

I Zform = 10, we integrate SFRD TO t and find a good agreement with the published SMDs. 

Conclusions. 
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1 . Introduction (FIR = bolometric IR). They find an increase of the fraction of 

hidden SFR from 56% locally to 84% at z = 1. The LDs show a 

One of the major objectives in astrophysics during the last 15 significant evolution as the FIR LD evolves faster than the FUV. 

years or so has been to follow the cosmic star formation rate den- Their ratio PfirIpfuv increases from ~ 4 ( A FUV ~ 1.3 mag) lo- 

sity (SFRD) at earlier and earlier epochs. But, whenever optical cally to ~ 15 (A FUV ~ 2.3 mag) at z = 1 . |Cucciati et al.(2012)1 

data are used, one must apply a dust correction to luminosity used the VIMOS- VLT Deep Survey to show, from the FUV 

densities (LDs) and a calibration into SFRDs (with their associ- only that the mean du st attenuation A FUV is in agreement with 

ated uncertainties) to obtain a relevant estimate. Knowing how |Takeuchi et al. (2005 )| over the range < z < 1. Then it remains 

the dust attenuation evolves in redshift is therefore mandatory if at the same level up to z ~ 2, and declines to ~ 1 mag at z ~ 4. 
one wishes to study the redshift evolution of the SFRD. 

For instance, Takeuchi et al.(2005) estimate the cosmic evo- In this letter, we use the FUV luminosity functions (LFs) 

lution of the SFRD from far-ultraviolet (FUV) and far-infrared published by |Cuccia ti et al.(2012) from the VLT along with the 
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FIR LFs fro m HerscAei-PACS+SPIR E datiQof a PACS selected 
sample from Gruppioni et al.(2013) to constrain the redshift evo- 
lution of log l0 (LFiR/Lfuv) (aka IRX) up to z ~ 4 for the first 
time directly using FIR data. With this information in hand, 
we can estimate the redshift evolution of PfirIpfuv as well as 
Ptot = Pfir+Pfuv- Finally, by integrating pror, we estimate the 
cosmic evolution of stellar mass density (SMD) with redshift. 

Throughout this paper we adopt a ACDM cosmology with 
(H ,Q. m ,Cl A ) = (70, 0.3, 0.7), where H is in kms^'Mpc^ 1 . All 
SFR and stellar masses presented assume, or have been con- 
verted to, a Salpeter IMF. 



2. Luminosity Functions 



Our analysis at z 



is based on the FUV LF from 



Wyder et al.(2005) and the FIR LF from Takeuchi et al.(2005) 



and for < z < 4 o n the FUV LF from |Cucciati et al.(2012)| 
and the FIR LF from |Gruppioni et al.(2013)| In the FIR and at 
z > 0, the sample is selected in the PACS bands but uses the 
full Herschel-PACS + Herschel-SPIRE SED data. The PACS 
selection means that we can miss sources towards the upper 
end of the redshift range. The LFs are evaluated from homo- 
geneous datasets in the FUV and the FIR. This minimizes bi- 
ases and keeps the same reference indicator throughout cosmic 
times with a simple well-defined and controlled selection func- 
tion. This is one of the strengths of this work. The FUV LFs 
are not corrected for dust attenuation. We define the LFs as a 
number density of galaxies with luminosity in logarithmic inter- 
vals, [logioL, logioL + dlogioL], where <I>(L) = dn/dlogioL and 
the luminosity is defined as L = vL v . FIR luminosities are de- 



fined as: <E>(L) 



<5*(^) 1 - Q; exp(- 



1 [log I0 d + L 



2,r- 



u 



)r). 



Observed uncertainties from |Cucciati et al.(2012)| and from 
Gruppioni et al. (20 13) are used whenever available. However, 



some of the Schechter parameters are fixed when the LFs are de- 
rived, namely a for the FUV LFs and a plus <x for the FIR LFs. 
Both in FUV and in FIR, we assume uncertainties of 10% up to 
z — 1, 20% up to z — 2 and 40% beyond for these fixed param- 
eters. This level of uncertainty is similar to previous works in 
FUV by e.g.|Oeschetal.(2010)[ |van den Burg et al. (2010) | and 
in FIR by Casey et al.(2012) We propagate uncertainties by sim- 
ulating 2000 realizations drawn from 1-cr Gaussian distributions 
for each parameter with known uncertainties and from a flat dis- 
tribution (i.e. equiprobability) for the fixed ones. We assume that 
all fixed values are equiprobable given the weak observational 
constraints. Finally, we interpolate the FUV and FIR Schechter 
parameters on the same redshift grid between z — and z — 3.6. 
Tab. [T] and Fig. [T] show the redshift variation of the LFs 
in FIR. The known difference in the FIR and FUV LFs 
|Takeuchi et al.(2005)| are clearly illustrated here: bright FIR 
galaxies are more numerous than bright FUV galaxies at 
logio(L[L Q ]) > 10. In FUV, except in the highest redshift bins, 
L* and <$>* remain approximately constant while the faint-end 
slope evolves. The FIR faint end slope is not observationally 
constrained at high z, and |Cucciati et al.(2012)| fix it to a - 1.2. 
However, L* and <D* are allowed to change with redshift. These 
different evolutions of the FUV and FIR LFs are reflected in 
Fig. [T] and explain the evolution of the cosmic SFRD and dust 
attenuation. 



1 From two Herschel Large Programmes: PACS Evolutionary Probe 
(Lutz et al.(201 1)) and the Herschel Multi-tiered Extragalactic Survey 



3. Dust Attenuation traced by the FIR to FUV LD 
ratio 

Fig. [2] presents the dust attenuation in the FUV vs. z and 
the ratio of the FIR to FUV LDs integrated in the range 
log w (L[L Q ]) = [7, 14] in the FUV (i.e. L'™ v = 1.65 x 
10~ 4 L* =3 , |Bouwens et al.(2009)) and [8, 14] in the FIR. The 
FUV dust attenuation is estimated from the IRX and converted 
to Afuv using Burgarella et al.(2005p] The redshift evolution 
of A FUV is in agreement with Cucciati et al.(2 012)[ Note that 
|Cucciati et al.(2012)| estimated Afuv through an analysis of in- 
dividual SEDs up to A bs - 2.2/mi (Ks-band). Fig. [2 suggests 
the presence of a local minimum at z ~ 2 that might be due to 
UV-faint galaxies (see Fig. 7 in |Cucciati et al.(2012) 1 that are 
responsible for a peak observed in the FUV LD and not ob- 
served in the FIR. Since the fields observed in FUV and in FIR 
are not the same ones, another origin might be found in cos- 
mic variance. The bottomline is that the existence of this trough 
in A F uv seems dubious. Finally, higher redshift A F uv from the 
UV slope, ft, suggest a continuous decline at least up to z = 6 
flBouwens et al.(2009)| >. 

We conclude that the cosmic dust attenuation Afuv reaches 
an absolute maximum at z ~ 1.2 followed by a global decline to 
z - 3.6 where it reaches about the same level measured at z - 0. 

The ft method is popular because estimates the total SFR 
from the FUV only. This is most useful at high redshifts 



where the samples are UV-selected (Burgarella etal. (2011) 
Bou wens et al.(2012)| |Heinisetal.(2013)| l. We propose to fol- 
low the redshift evolution of the cosmic volume-averaged points 
in the IRX - ft diagram (Fig. [3T> to constrain models. However, 
we must caution that the values plotted in Fig. [3] cannot be di- 
rectly compared to galaxies. The x-axis is calculated from the av- 
eraged rest-frame FUV - near-UV colors ( |Cortese et al.(2006)) . 
Horizontal error bars indicate the dispersion of the FUV slope. 
The IRX is estimated from LFs and is therefore volume- 
corrected. Vertical error bars are uncertainties. This IRX -ft plot 
can be interpreted as the location of a comoving volume as a 
function of redshift. From z ~ 1 to z ~ 4, the points evolve 
downwards pa rallel to the original |Meurer e t al.(1999) law and 
the update by |Takeuchi et al.(2012)| 

4. The total FUV+FIR Star Formation Density and 
Stellar Mass Density 



The calibration from 



LD to SFRD is 
n the FUV and also 



problematic 
in the FIR 



( Kobayash ietal.(2013) i_ 

(Kennicutt(1998), Schaerer et al.(2013)). In agreement with 



|Casey et al.(2012) 



Kennicutt(1998) calibrations and 



assume a Salpeter initial mass function (IMF) to allow a better 
comparison with other published SFRDs. Note that the AGN 
contribution to the FIR LDs has been estimated in FIR using a 
SED fitting and subtracted off i.e. the presented FIR LD is due 
to star formation only. 

Fig. [4] suggests a flattening of the total SFRD up to 
3 (|Chary & Elbaz(2001)| [Perez-Gonzalez et al.(2005")l 



Le Floc'h et al.(2005)| |Franceschini et al.(2010)) where the UV 
data favour a peak followed by a decrease. Note that we could 
not rule out a small increase or decrease within the uncertain- 
ties. The plateau up to z ~ 2.5 is seen by |Rodighiero et al.(2010)| 
and Magnelli et al.(2011) while the decrease of Psfr at z > 



2 The conversion from IRX to Afuv from Burgarella et al.(2005) is 



(HerMES, |Qliver et al.(20T2J} 



valid at log 10 (L PIR /L PUV )>-1.2: A F 



-0.028 [log 10 L FIR /L FUV Y + 



0.392 [log 10 L FIR /L FUV Y + 1.094 [log 10 L FIR /L FUV ] + 0.546 
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Table 1. Schechter parameter for FUV and FIR luminosity functions. 



Redshift range 


M* or L* fl 


o* 


a 


CT b 


FUV luminosity functions 


0.0<z<0.1 c 


-18.04 ±0.11 


-2.370 ± 0.06 


-1.22 ±0.07 





0.05 < z < 02 d 


-18.12 + 0.00 


-2.155 + 0.03 


-1.05 ±0.04 





0.2 <z< 0.4 rf 


-18.3 + 0.20 


-2.161 ±0.06 


-1.17 ±0.05 





0.4 < z < 0.6 rf 


-18.4 + 0.10 


-2.180 + 0.06 


-1.07 ±0.07 





0.6 < z < 0.8 d 


-18.3 + 0.10 


-2.021 ± 0.05 


-0.90 ± 0.08 





0.8 <z< 1.0 d 


-18.7 + 0.10 


-2.045 ± 0.05 


-0.85 ±0.10 





1.0 <z< I2 d 


-19.0 ± 0.20 


-2.129 + 0.07 


-0.91 ±0.16 





1.2 <z< l.l d 


-19.6 ± 0.20 


-2.387 + 0.10 


-1.09 ±0.23 





1.7 < z < 2.5 d 


-20.4 + 0.10 


-2.472 ± 0.03 


-1.30 ±-0.26 





2.5 < z < 3.5 d 


-21.4 + 0.10 


-3.066 ± 0.03 


-1.50 ±-0.60 





3.5 < z < 4.5 d 


-22.2 ± 0.20 


-3.959 ± 0.04 


-1.73 ±-0.69 





FIR luminosity functions 


z = e 


9.25 ± 0.00 


-2.051 ±0.00 


1.23 ±0.00 


0.72 ± 0.00 


0.0 < z < 0.3 / 


10.12 + 0.16 


-2.29 ± 0.06 


1.15 ±0.05 


0.52 ± 0.05 


0.3 < z < 0.45' 


10.41+ 0.03 


-2.31 ± 0.03 


1.2 ±-0.12 


0.5 ± -0.05 


0.45 < z < 0.6-f 


10.55+ 0.03 


-2.35 ± 0.05 


1.2 ±-0.12 


0.5 ± -0.05 


0.6 < z < 0.8 / 


10.71 ± 0.03 


-2.35 ± 0.06 


1.2 ±-0.12 


0.5 ± -0.05 


0.8 < z < l.Of 


10.97 ± 0.04 


-2.40 ± 0.05 


1.2 ±-0.12 


0.5 ± -0.05 


1.0 < z < 1.2 / 


11.13 + 0.04 


-2.40 ± 0.05 


1.2 ±-0.24 


0.5 ±-0.10 


1.2 < z < 1.7 / 


11.37 + 0.03 


-2.70 ± 0.04 


1.2 ±-0.24 


0.5 ±-0.10 


1.7 < z < 2.0 / 


11.50 + 0.03 


-2.85 ± 0.03 


1.2 ±-0.24 


0.5 ±-0.10 


2.0 < z < 2.5^ 


11.60 + 0.03 


-3.01 ±0.11 


1.2 ±-0.48 


0.5 ± -0.20 


2.5 < z < 3.0^ 


11.92 + 0.08 


-3.27 ±0.18 


1.2 ±-0.48 


0.5 ± -0.20 


3.0 < z < 4.2/ 


11.90 + 0.16 


-3.74 ±0.12 


1.2 ± -0.48 


0.5 ± -0.20 


Notes. The top panel lists the FUV LFs and the bottom panel lists the FIR LFs. For all the p< 
assumed 20% or error. (a) L* [L Q ] for FIR LFs or M* [AB mag] for FUV LFs (,,) a only needed fc 


irameters with an uncert 
>r FIR LFs <c) from|Wyde 

(222,0-! = 0.5135 +0.( 
94 ± 0.0222, o- 2 = 0.8 
igher redshifts, we mac 


ainty set to 0.00. we 
r et al.(2005)| w from 


|Takeuchi et al.(2005)| (eJ |Cucciati et al.(2012)|^|Grup 


pioni et al.(2013)| 

nodel based on witha; = 0. 1261+0. C 
3d from LBGs. 0.0959 and a 2 = 0.22 
nt with that of 2.7151 ±0.0839. Ath 
ange. However, are explained below. 




2.5 is predicted by Betherr 
evolution of mass function 
All in all, our total SFRD 


nin et al.(2012a)|s r 
and sSFR estimat 
is in fair agreeme 

i the same redshift t 


(704, z\ = 1.1390± 
160 ± 0.0964, z 2 = 
Je assumptions that 


Hopkins & Beacom(2006)|ij 





discrepancies exist: our total SFRD is lower at z < 1 and is 
only marginally consistent but lower at z > 3. Also, note that 
PACS data are less sensitive at higher than at lower redshift since 
the rest-frame wavelength moves into the mid-IR. The prelimi- 
nary FIR SFRD from |Vaccari et al.(2013)1 (Hersriiei/SPIRE se- 
lection) is in excellent agreement over the < z < 2 range 
but is slightly higher than that derived from PACS at z > 3. 
However, this is only a ~ 2<x difference. Barger et al.(2012) pub- 
lished a FIR SFRD based on SCUBA-2 data that is also in agree- 
ment with ours at 2 < z < 4. We first tried to fit SFRDjot 



with a one-peak analytical function (Hopkins & Beacom(2006) 



Behroozi, Wechsler & Conroy(2012) i but the results are not sat- 



isfactory. So, we combined two Gaussians: 



-I.-.-2'- 



a\e 



+ aie 



The cosmic SFRD presents a (weak) maximum at z ~ 
2.5 - 3.0 (i.e. between 2.6 - 2.1 Gyrs resp.) while the dust at- 
tenuation presents a maximum at z ~ 1.2 (i.e. 5 Gyrs). We 
have tried to lock the faint end slope of the UV LF -1.2, to see 
how far out in redshift the peak of obscuration could potentially 
move and we do not detect any change, suggesting this effect 
is solid. We have no definite explanation for this delay of ~ 2.7 
Gyr. Type II supernovae start producing dust earlier than AGB 
stars ( Valiante et al.(2009)) but the difference in timescales is too 
short and only on the order of a few 10 Myr for the onset of dust 
formation. Dust grain destruction in the ISM might play a role 
( |Dwek & Cherchneff(201 1)) but the efficiency of destruction is 
not well known and depends on the star formation history. These 
dust-related origins for the delayed maximum are not likely. The 
best explanation might be that this delay is related to a global 
move of galaxies in the [log I0 (L f// {/L F£/ y) vs. log 10 (L F /« + 
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Fig. 1. Redshift evolution of the FIR from Grup pioni et al.(2013)| (red) and FUV from Cucc iati et al.(2012)] (blue) LFs. Note that the 
FUV LFs are uncorrected for dust attenuation. The LFs at every other redshift are plotted boldly. The others are fainter to lighten 
the figure. The LFs are plotted within the limits of integration. 



Lfuv)] diagram. JBuat et al.(2009) showed that galaxies evolve 
in redshift from z = to z = 2 in this diagram, with high red- 
shift sources having lower IRX at given total luminosities. This 
change is likely related to systematic changes of the FIR SEDs 
themselves ( [Elbaz et al.(2011)| |Nordon et al.(2012)| >. This sug- 
gests that the shift might be caused by the relative importance of 
more luminous galaxies (log l0 Lfuv[L Q ] > 10) in the FUV as z 
evolves. 

By integrating the SFRD, we can estimate the stellar-mass 
density (SMD) (Fig. B] and Tab. |2l). To do so, we stress that 
we set the mass fraction of a generation of stars that is re- 
turned to the interstellar medium to a fixed value R = 0.3 
(Fraternali & Tomassetti(2012)). We also have to assume a star 
formation history from z - 3.6 up to the galaxy formation set 
at Zform - 10. Option 1 is a linear extrapolation while option 
2 corresponds to a rising exponential e'^ T with t = 0.42 as 
in Papjovich et al.(2011) that joins the observationally-deduced 



'apo\ 



SFR 

Superimposed in Fi g. [5] are recent S MDs (converted to 
Salpeter IMF if needed). Stark et al.(2013) account for the neb- 
ular emission lines contribution to the broad-band fluxes used 



Selecting option 1 or 2 does not impact on Fig. 4 



to infer stellar masses (Ono et al.(2010), de Barros et al.(2012)). 
The trend from Labbe et al.(2 010) lies above our points. The oth- 
ers SMDs are in agreement within the uncertainties at 0.6 < z < 
3.6 for the two above options. Wil kins et al.(2008J1 compile mea- 
surements of the SMD from the literature and provide a best fit 
parametric law p*(z) = ae~ hz ' where a - 0.0023, b - 0.68 and 
c — 1.2. We also overplot it in Fig.B] This curve slightly under- 
estimates our SMF at very low redshifts but at higher redshifts, 
it follows the points derived from our data and our assumptions. 

We reach a fair agreement, especially at 0.6 < z < 3.6. 
The discrepancy previously observed is reduced here but still 
marginally consistent at very low redshifts. As shown in Fig. |4| 
our total SFRD generally lies below Hopkins & Beacom(2006) 
suggesting that this previous evaluation of dust attenuation 
might be over-estimated in this redshift range. Making use 
of FIR data allows us to reach a better agreement. Note that 
Hopkins & Beacom(2006) did not directly use MIR-based data 
to estimate their best-fitting parametric curve but only to correct 
the obscuration of the FUV data. 

Fig. [6] shows the variation of the cosmic specific star for- 
mation rate (sSFR = SFRD/SMD) as a function of redshift. 
Previous observational results often suggest a flat evolution of 
the sSFR at z > 2 ( |Gonzalez et al.(2012)irBouwens et al.(2012)| 



Burgarella, et al.: < z < 4 redshift evolution of the UV+IR SFRD and dust attenuation 




A FUV this work 
best fit V:^> 

l+(z/c) 

Cucciati et al. 2012 
Takeuchi et al. 2005 
Bouwens et al. 2009 



2.5 



2.0 



> 

fa 

1.5 < 



1.0 



0.5 



redshift 

Fig. 2. Left axis: ratio of the FIR to FUV LDs (IRX). Right axis: FUV dust attenuation (A/rj/y). The red dotted line with red diamonds 
is from Takeuchi, Buat & Burgarella (2005). The green filled area and green dots are the associated uncertainties estimated through 
bootstraping with 2 000 drawings. Black dots are the values directly computed from the LFs. At z — 3.6, Apuv reaches about the 
same value as a tz — 0. [T akeuchi et al.(2005) (red diamonds) used an approach identical to ours while a SED analysis (no FIR data) 
is performed in Cucciati et al.(2012) (blue boxes). Bouwens et al.(2009) are estimates based on the UV slope /?. The limiting FUV 
luminosity is 10 7 L o or 1.65 X 10~ 4 Z/ 3 . The best fit is given by A FUV (z) = -^^ with a = 1.20, b = 1.50, c = 1.77, and d = 2.19. 



|Schaerer etal.(2013)|l while most the oretical models predict 



a continuous rise ( Bouche et al.(2010j] |Weinmann et al.(201 1)| 
|Daveet al.(2011)| i. No matter what hypothesis selected to ex- 
trapolate the observed SFRD beyond z > 3.6 to Zform = 10, 
our sSFRs remain consistent with an increase at low redshifts. 
The influence of the SFRD assumed at z > 3.6 is not notice- 
able within the uncertainties at z < 3.6. A comparison with the 
sSFR of galaxies f rom jNoeske et al.(2007)l |Daddi et al.(2007)| 
Wuyts et al.(2011) and |Bouwens et al.(2012) in Fig. 16] suggests 



that the most massive galaxies {log^M* — 10 -10^M G ]) are 
in agreement with our sSFRs (we corrected the SMDs to R — 0.3 
and we applied a correction to the calibration to SFR if neces- 
sary). Finally, we note that at higher redshifts, option 1 keeps on 
rising while option 2 shows a flattening at 4 < z < 5 followed 
by an increase at z > 5. We stress, though, that by assuming an 
exponential rise above z > 4 with a value of the time constant 



t = 420 Myrs as in Papovich et al.(2011) and the formation of 
galaxies at Zf orm =10 implies that the observed plateau at 4 < z < 
5 must be temporary. Changing t and/or Zf orm to larger redshifts 



would shift the increase in sSFR to earlier times. Theoretically, 
fixing Zf orm = oo would translate into a flat sSFR. 



5. Discussion and Conclusions 

The variation of the cosmic dust attenuation with redshift as es- 
timated from the IR to FUV luminosity ratio suggests the pres- 
ence of a peak in the dust attenuation at z ~ 1 .2 followed by a 
decline up to z — 3.6. This result confirms |Cucciati et al. (20 12)1 
dust attenuation estimated from SED fitting without FIR data. 
Moreover, the redshift evolution of the volume-corrected IRX -ft 
points globally follows the IRX - /3 law from z = 0.4 to z — 3.6. 
The total (FUV+FIR) cosmic SFRD starts declining above 
z — 3 - 4 and reaches the same level at z ~ 5 - 6 as is measured 
locally if we assume no variations in this trend. At z - 3 - 4, the 
decrease observed in the SFRD is not unexpected: most high red- 
shift studies clearly suggest such a trend through rest-frame UV 
observations (Hopkins & Beacom(2006) l Bouwe ns et al.(201 l)| l 
and predicted by Bethermin et al. (2012a) Backwards in time 
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Fig. 3. Dust attenuation vs. redshift. The boxes are color-coded according to redshift. Note that the x-axis bars correspond to 
the dispersion in UV slope while the y-axis are evaluated from the uncertainties in the LFs. The black dots are the original 
data points from |Meurer et al.(1999)| and the black curve is Meurer et al.(1999 j]s law. T he black dashed line and grey boxes 
are the update (using the same apertures in FUV and in FIR) from Takeuchi et al.(2012)l Strictly speaking, our points and 
those from Meurer et al.(1999) are not comparable because we use volume-corrected LFs and not individual galaxies as done 
in Meurer et al.(1999) In the diagram, they show an almost continuous decrease with increasing redshift and lie in between the 
Meurer et al.(1999) and Tak euchi et al.(2012)| laws. It is important to stress that the dust attenuation A FUV are estimated from the 
IRX and not from the U V slope j3. 



from today, this decline is preceded by a rise from z ~ to a 
break at z ~ 1-2, followed by a plateau up to z — 3 - 4. 

If we compare Fig. [2]to Fig. HI the peak of the dust attenua- 
tion is delayed with respect to the plateau of the total SFRD. It 
is also coeval with the final decrease at z = 1 — 1.5, the peak 
of the less luminous AGNs (Hopkins(2007)). A similar peak 
seems to appear at z - 1-2 (depending on the stellar mass) 
of the co smological merger rate ([C onselice et al.(2008) ) and of 
the CIB ( |Bethermin et al.(2012b)) . Are all these effects related 
to the same physical phenomena and what are their characteris- 
tic timescales? To better understand the delay, it is necessary to 
perform an analysis via models that are fed with data of the gas 
content and the evolution of metallicity. 

Using the observed cosmic SFRD along with the assumption 
of an exponential rise from Z/ orm = 10 to z — 3.6, we are able to 
recover the SMD evaluated from galaxy surveys. With the same 
assumption, we predict a flattening of the sSFR at 3 < z < 5 
followed by a new steepening at z > 5. 



Fig. [2] and Fig.ffltaken together at face value would suggest 
that the universe's dusty era (meaning dust attenuation larger 
than in the local universe) started at z = 3 - 4 simultaneously 
with the rise of a universe-wide star-formation event. 

Fig. |4] Fig. [5] and Fig. [6] allow us to follow the 
SFRD, the SMD, and the sSFR over most of the Hubble 
time in a consistent way. However, large uncertainties pre- 
vent us from closing the case. Additionaly, it remains quite 
puzzling that GRB-based analyses suggest a much shal- 
lower decrease ( |Kistler et al.(2009)| |Robertson & Ellis(2Q 1 2)] > 
than Lyman break galaxies. The statistical significance of 
these results is still debated because of the low number 
of objects at high redshifts and a possible modification of 
the IMF ( |Dwek & Cherchneff(201 1)1 |Haywardetal.(2013)) . 
Another possibility is that GRBs might still be biased toward 
certain types of SFGs, even though this bias may be less than 
thought a few years ago. 
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Fig. 4. SFRD densities in the FUV (blue), in the FIR (red) and in total (i.e. FUV + FIR) in green (other colors are due to overlaps of 
the previous colors). The lines are the mean values while the lighter colors shows the uncertainties evaluated from the 2 000 runs as 
in Fig. [2 After the initial increase of the total SFRD from z — to z ~ 1.2, it remains flat or slightly increases/decreases up to z 



2.5-3.0 followed by a decrease. Globally and over < z < 3.6, the total average SFRD is slightly below Hopkins & Beacom(2006) 
and in agreement with Behroozi, Wechsler & Conroy(2012) up to z ~ 2.. SFRD from Barger et al.(2012)Jand preliminary results 
from Herschel/S PIRE estimated by |Vaccari et al.(2013)| are in agreement with these trends. Symbols and lines are explained in the 
plot. 
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Fig. 5. Evolution of the cosmic SMD vs. redshift. The lines with blue and red boxes are the mean values evaluated by integra- 
tion while the lighter colors show the uncertainties from the 2 000 runs as in previous figures. We overplot the trends found by 
|Labbeetal.(2010)| at z-1 - 8, |Gonzalez et al.(2011)|and|Stark et al.(2013)l at 3 < z < 7. We also plot the line corresponding to the 
compilation of published measurements by Wilkins et al.(2008) Within the uncertainties, we find a good agreement of the SMD 
integrated from the SFRD with all other SMDs based on galaxy surveys from z - 0.6 to z = 5. The black- limited area represents 
the compilation of results by Santini et al.(2012) Scaled to the same cosmology and IMF (Salpeter). Using Gonzal ez et al.(201 1)| 
we integrate down to luminosities equivalent to M ~ 2500M o . 
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10.5 - 11.0 galaxies. At z > 4 the extrapolations are in agreement with log 10 M* ~ 9.5 as evaluated 



10 



